Low-Dimensional Long-Range Topological Charge Structure in the QCD Vacuum 



O 

o 

(N 



I. Horvath", S.J. Dong", T. Draper*^, F.X. Lee'''^ K.F. Liu°, N. Mathur", H.B. Thacker'' and J.B. Zhang'^ 
'^Department of Physics, University of Kentucky, Lexington, KY 40506 
Center for Nuclear Studies and Department of Physics, George Washington University, Washington, DC 20052 
"Jefferson Lab, 12000 Jefferson Avenue, Newport News, VA 23606 
'^Department of Physics, University of Virginia, Charlottesville, VA 22901 
" CSSM and Department of Physics, University of Adelaide, Adelaide, SA 5005, Australia 

(Dated: September 15, 2003) 

While sign-coherent 4-dimensional structures cannot dominate topological charge fluctuations in 
the QCD vacuum at all scales due to reflection positivity, it is possible that enhanced coherence 
exists over extended space-time regions of lower dimension. Using the overlap Dirac operator to 
calculate topological charge density, we present evidence for such structure in pure- glue SU(3) lattice 
gauge theory. It is found that a typical equilibrium configuration is dominated by two oppositely- 
charged sign-coherent connected structures ( "sheets") covering about 80% of space-time. Each sheet 
is built from elementary 3-d cubes connected through 2-d faces, and approximates a low-dimensional 
curved manifold (or possibly a fractal structure) embedded in the 4-d space. At the heart of the 
sheet is a "skeleton" formed by about 18% of the most intense space-time points organized into a 
global long-range structure, involving connected parts spreading over maximal possible distances. 
We flnd that the skeleton is locally 1-dimensional and propose that its geometrical properties might 
be relevant for understanding the possible role of topological charge fluctuations in the physics of 
chiral symmetry breaking. 
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PACS numbers: ll.15.Ha, 11.30.Rd 

One of the important goals in hadronic physics is a 
detailed understanding of vacuum structure in QCD. 
Local patterns in topological charge fluctuations repre- 
sent an important aspect of this structure. Indeed, phe- 
nomena such as large 77' mass, ^-dependence, and possi- 
bly spontaneous chiral symmetry breaking (SxSB) are 
directly related to vacuum fluctuations of topological 
charge. While lattice gauge theory provides a framework 
for first-principles non-perturbative studies of topologi- 
cal charge fluctuations, the extraction of meaningful and 
unbiased structural information has long been a difhcult 
problem. However, the recent development of fermionic 
methods has improved this situation considerably and it 
is no longer necessary to rely on subjective procedures to 
smooth out the rough short-distance behavior of gauge 
fields. One way to proceed is to study the structure of 
low-lying Dirac eigenmodes and to infer the properties 
of underlying topological charge fluctuations indirectly. 
This was advocated in P], developed in 0, and used in 
many recent studies (see e.g. 

The advances in implementing lattice chiral symme- 
try (for reviews see e.g. ^) have also made it possible 
to pursue a direct approach to this problem. Indeed, 
with any 75-IIermitian lattice Dirac operator D satisfy- 
ing the Ginsparg- Wilson relation {D,"f^} = D"f^D (e.g. 
the overlap operator 0)jOne can associate a topological 
charge density operator [7| 
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The global charge associated with this lattice density is 
strictly stable with respect to generic local variations of 
the gauge field 0], thus fully respecting the topological 
nature of this quantity. In addition, the above qx is spe- 



cial in several ways. For example, the index theorem 
is exactly satisfied with respect to the associated Dirac 
operator 1^ , and its renormalization properties are analo- 
gous to those in the continuum 9]. An important virtue 
of qx is that it can be naturally eigenmode-expanded. 
While qx contains all fluctuations up to the lattice cutoff, 
one can use the eigenmode expansion up to scale A to de- 
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fine an effective density [iflj , q, 

where = V'^^TsV'^ is the local chirality of the mode 
with eigenvalue A. The ultraviolet fluctuations are natu- 
rally filtered out in q'J^^ {eigenmode filtering 2]), and this 
density can be used to study the structure of topological 
charge fluctuations at arbitrary low-energy scale A. By 
studying the effective densities, it was demonstrated [I^ 
that topological charge at low energy is not locally con- 
centrated in sign-coherent unit-quantized lumps (e.g. in- 
stantons). This was first predicted by Witten [ll| and 
emphasized recently in , where the first lattice evidence 
was presented. 

In this work, we initiate the investigation of possible 
non-trivial structure (order) in the full topological charge 
density Eq. for typical configurations contributing to 
the pure glue QCD path integral. This is a qualitatively 
new step with some aspects worth emphasizing, (a) We 
evaluate the local operator qx at every space-time point, 
thus producing gauge invariant "configurations" of topo- 
logical charge density. No processing of the underlying 
gauge configurations is involved and no bias is imposed, 
(b) Assuming that the well-defined space-time structure 
in qx exists, it has to be viewed as fundamental (rather 
than low-energy) since, unlike the case of effective den- 
sities q^\ the use of a local operator does not intro- 
duce a new scale apart from the lattice cutoff already 
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in place. Consequently, such structure could have both 
short and long-distance manifestations, (c) A fundamen- 
tal structure of this type has not been observed before 
in an unbiased setting. Indeed, using the conventional 
operators, the resulting space-time distributions of topo- 
logical charge can hardly be distinguished from complete 
disorder. The common interpretation of this is that the 
physically relevant fluctuations might be obscured by the 
structureless ultraviolet noise arising due to entropy con- 
siderations (the problem of ultraviolet dominance). 

In what follows, we will present evidence for the ex- 
istence of fundamental structure in topological charge 
density, and demonstrate some of its basic geometric fea- 
tures. Needless to say, we do this with the aim that 
some of these properties might help to advance our un- 
derstanding of the QCD vacuum, and hence of the basic 
features of hadronic physics. With regard to the problem 
of ultraviolet dominance, we were motivated by the fact 
that, apart from beautiful properties mentioned above, 
the operator Qx of Eq. differs from standard lattice 
discretizations in a very qualitative manner. In particu- 
lar, it is well-known that lattice chiral symmetry implies 
non-ultralocality of lattice Dirac kernel D, i.e. non-zero 
coupling among fermionic variables at arbitrarily large 
distances Similarly, it is believed that any D^y in- 
volves gauge paths (from x to y) running arbitrarily far 
from X and y. As a result, qx will receive (small) con- 
tributions from arbitrarily extended gauge loops {chiral 
smoothing I3| ) and is expected to be much less sensitive 
to the ultraviolet noise than typical ultralocal operators. 
We find that chiral smoothing indeed provides enough 
suppression of ultraviolet noise for the underlying coher- 
ent structure to be revealed. 

As a guide to identifying this structure, we note the 
fact that ( q{x)q{0) ) < , |a;| > , in the continuum [l3l |. 
This implies that the topological charge cannot be pre- 
dominantly concentrated in 4-dimensional sign-coherent 
structures of finite physical size However, one can- 
not conclude that the local behavior of topological charge 
density is strictly disordered! Indeed, the negativity of 
the correlator does not exclude enhanced sign-coherence 
present on lower- dimensional subsets of 4-d space-time. 
In that case the structure could respect the negativity of 
the correlator by appropriate embedding of subsets with 
alternating sign into space-time. Consequently, we fo- 
cus our search on the level of sign-coherence, and how 
it changes when concentrating on space-time subregions 
containing strong fields. In what follows we will refer to 
sign-coherence in qx simply as coherence. 

For the numerical study we use qx associated with the 
overlap Dirac operator 6] , whose suitability for this pur- 
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TABLE I: Ensembles of Wilson gauge configurations. The 
global topological charges are (3,0,-2,0,2,-1,-1,3) for £\ and 
(2,1) for £2. 



pose was confirmed in Ref. m . To sample the QCD vac- 
uum, we work with Wilson gauge backgrounds as summa- 
rized in Tabled Configurations are separated by 20,000 
sweeps and the Sommer parameter serves to set the scale 
using the interpolation formula of Ref. Since the 

computation is extremely time-consuming, the ensembles 
are small. However, the attempts to understand QCD 
vacuum structure in the path integral framework implic- 
itly rely on the expectation that relevant individual con- 
figurations exhibit common qualitative properties that 
are "typical" of the QCD vacuum, so that interesting re- 
sults should be obtainable from even a few configurations. 
In our case, we found that for the aspects discussed be- 
low, the behavior is remarkably stable from configuration 
to configuration (both qualitatively and quantitatively). 
In fact, one could infer the main results presented here 
by studying just a single configuration from each ensem- 
ble. We are thus confident that our conclusions are not 
affected by low statistics. Some preliminary results re- 
lated to this work were presented in Ref. ■ Details of 
the fermion implementation can be found in Ref. 0- 

1. We start by verifying that qx is indeed not dom- 
inated by coherent 4-d structures, as argued above. 
Regardless of their shapes, the possible candidates for 
such structures must be built of elementary lattice 4-d 
cubes 17] that are themselves coherent (their sites have 
the same sign of qx). In fact, we can define the coherent 
regions on the lattice by finding all coherent 4-d cubes Cx 
(labeled by origin x) and identifying maximal connected 
regions built from such cubes. We will say that lattice 
region TZ = {Cx,]i = 1, . . . , iV } is connected, if for arbi- 
trary Cxi , Cxj € TZ there is a sequence Cy^. , A: = 1, . . . , n , 
such that (a) Cy^. £ TZ, (b) yi — Xi, yn — Xj, and (c) Cy^. 
and Cy^j^-^ share a common face (3-d cube). 

We have determined all such regions ( "structures") TZk 
present in a given configuration, and ordered them by the 
number of lattice sites in the structure, Nk- If there are 
4-d structures of finite physical size dominating in the 
continuum limit, then at least the largest of these co- 
herent regions should exhibit scaling behavior when the 
lattice spacing is changed. Since the ensembles £1, £2 
have the same physical volumes, the ratios N^/V {V is 
the number of lattice points) for typical largest struc- 
tures should not depend on the lattice spacing. In Fig.Q] 
(top) we show these ratios for a configuration from both 
ensembles. Rather than being similar, the ratios drop 
dramatically at the smaller lattice spacing. To see this 
on average, we compute the mean value of the ratio in- 
cluding the 5 largest structures per configuration. Fig. ^ 
reveals that the ratio not only drops, but it drops by a 
factor larger than the ratio of lattice volumes. This is 
confirmed by the fact that the average number of 4-d 
cubes per large structure is about 35.3 at a = 0.110, and 
is about 29.4 at a = 0.082. The volume of structures thus 
decreases even in lattice units. These considerations not 
only confirm that 4-d structures do not dominate topo- 
logical charge fiuctuations, but also suggest that such 
structures might not occur in the continuum at all. In- 
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FIG. 1; Top: The fraction of space-time points occupied by the largest 20 structures buih of 4-d cubes on a typical configuration 
from ensemble £i (left) and £2 (middle). The average fraction from the 5 largest structures is also shown (right). Bottom: The 
same as top but for structures built of 3-d cubes. There are error bars on the rhs plots. 



deed, in physical units, structures TZk observed on the 
lattice appear to shrink to mere points in the continuum 
hmit. 

2. We now concentrate on the possibility that qx 
exhibits coherence on lower-dimensional subsets of 4-d 
space-time. As a first step, we search for maximal coher- 
ent regions built from elementary 3-d cubes connected 
through 2-d faces. Such regions define connected 3-d lat- 
tice hypersurfaces embedded in the 4-d space-time. We 
find that all of the configurations studied (see Table P) 
exhibit a remarkably similar coherent behavior. The typ- 
ical situation for a configuration from ensembles £1, £2 
is shown in Fig. ^ (bottom). Contrary to the case of 4-d 
cubes, there are only two (oppositely charged) coherent 
regions dominating the behavior of q^- Together, these 
folded lattice "sheets" cover about 80% of sites and prac- 
tically fill the space-time. This situation is insensitive to 
the change of the lattice spacing and our data indicate 
that it will survive the continuum limit (bottom right of 
Fig.HJ). 

We wish to emphasize two points. First, we have tested 
whether the existence of the 2-sheet structure reflects a 
specific order present in q^, or if it could occur in a ran- 
dom situation as well. To do that, we have performed 
a random permutation of sites p{x) on our lattices and 
looked for the coherent behavior in qp(^x)- The typical re- 
sult of such a test is shown in Fig. [21 The 2-sheet struc- 
ture disappears and the fractions involved are almost in- 
visible on the scale of the original plot. Secondly, the 
existence of (almost) space-filling 3-d lattice sheets is not 



sufficient to conclude that the local dimension of maxi- 
mal coherent regions in the continuum limit is three. The 
determination of this dimension requires a scaling anal- 
ysis which is in preparation. While possibly less than 
three, the dimension is at least one (see part 4 below). 
Note that this issue is not addressed by the lower right- 
hand plot of Fig. n There the implied scaling relates to 
the fact that the sheets occupy a macroscopic fraction 
of space-time and hence involve an infinitely large lower- 
dimensional manifold in the continuum limit (in the 4-d 
volume of finite physical size). 
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FIG. 2: A typical coherent structure before and after the ran- 
dom permutation of sites. On the common scale, the fractions 
on the r.h.s. plot are barely visible. Configuration from £2- 

3. Since the 2-sheet structure almost fills the space- 
time, it is expected that for a generic point a: on a given 
sheet, there are points y on the same sheet separated from 
X by maximal space-time distances possible. Because the 
sheet is path-connected, such points can be connected by 
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paths within the sheet. Thus, from the geometric point of 
view there is no specific long-distance scale (other than 
infrared cutoff) associated with global behavior of the 
sheet. We will refer to such a situation occuring in the co- 
herent structure as the super-long-distance property |l8j| . 

The super-long-distance property might be relevant to 
the physics of the QCD vacuum, e.g. to Goldstone boson 
propagation, if it remains true also if one concentrates 
on regions with largest topological charge density. In- 
deed, the physics of the vacuum is expected to be mostly 
driven by fluctuations with intense fields. For a given 
configuration, let iS-'^ (0 < / < 1) be the subset of the 
lattice containing the fV most intense sites ranked 
by 1 9a; I- Within S-^ , there are maximal coherent regions 
Wl'^ C S^, k ^ 1,2, ...,iV'*'/, built from d-dimensional 
cubes connected through (d— l)-dimensional faces. For 



G one can find a maximal Euclidean distance r 



f 



reachable from x by traveling on a path within a single 
structure T?.^'-'^ 9 x. Note that r;J ^ = for points that do 
not belong to any structure. Considering a fixed small 
fraction /o < /, we calculate the average of ^ over 
X £ S-^" , i.e. r^°'-^ = {r^^)gfo- This provides us with 
the measure (scale) of how extended the coherent struc- 
ture typically is. In Fig. O we show the /-dependence 
of r'^'^^'^ measured in units of the largest Euclidean dis- 
tance {Li are the sizes of the lattice box). 
The average is taken over the ensemble £i and the ob- 
served behavior is typical (and stable) for /o of a few 
percent. Also shown is the result after performing the 
random permutation of sites for every configuration. 

Note that in the case of 4-d cubes the structure never 
becomes super-long-distance. In fact, we expect that 
the maximal distance will approach zero in the contin- 
uum limit at arbitrary /. For d = 3, 2, 1, we observe 
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FIG. 3: The /-dependence of r^J ''^'^ for ensemble £i (up- 
per curves of each plot) and after random permutation of 
sites (lower curves of each plot). Results using 1,2,3 and 4- 
dimensional hypercubes are shown. 
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FIG. 4: Lattice coherence dimension. Explanation in text. 

a well-defined transition to a super-long-distance state 
which is certainly present at / « 0.8, 0.4, 0.2 respec- 
tively {r^'^^''^ ~ 0.9 at those fractions). These values 
are insensitive to the change of the lattice spacing. We 
thus see that after removing the low-intensity points, the 
super- long-distance structure remains at the core if lower- 
dimensional elementary cubes are used. The minimal 
structure with this property is exposed with 1-d cubes 
(links) and we will refer to it as the skeleton. In our nu- 
merical experiments, the maximal space-time distances 
start occuring already at about / = 0.16 while the aver- 
age saturates at about / = 0.20. We thus take / = 0.18 
as our initial reference value for the skeleton. As further 
evidence that the observed structure carries a significant 
amount of order, we note that at / = 0.18 (with the 
skeleton already formed) one is still confined to average 
maximal distances of a single lattice spacing when ran- 
dom permutation of sites is performed. 

4. Even without an extensive scaling analysis, it is pos- 
sible to obtain useful insight into the question of the ac- 
tual space-time dimension involved in the skeleton. More 
precisely, we will discuss the local dimension. To appre- 
ciate the difference, consider e.g. a network built from 
1-d lines, but organized globally into a 2-d structure 
("mesh") with a typical physical distance I among the 
nodes. Locally, such a manifold is 1-dimensional, but on 
the scales larger than I it will appear 2-dimensional. We 
will consider the global structure from this point of view 
elsewhere. 

Our tool here will be the lattice local dimension d^ iS) 
which we define for an arbitrary subset S of the lattice. 
For a; G 5 let d^^ be the maximal dimension d of the ele- 
mentary cube Cy such that x £ Cy and Cy C <S, i.e. the 
maximal dimension of the cube within S containing x. 
Then d^{S) = {d^)s is the average of this local dimen- 
sion over S. While d^{S) is strictly a lattice notion, it 
can still be quite useful. In particular, if there is an un- 
derlying set S'^ in the continuum characterized by a single 
local dimension d{S'^) < 4, then the continuum limit of 
d^{S) typically provides an upper bound for d{S'^). For 
example, consider a line of non-zero physical length in 
the continuum. At sufficiently small lattice spacing this 
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FIG. 5; Generic behavior of on a 2-d section of space-time 
(configuration from £i). 

might appear as the connected structure built predom- 
inantly of bonds, in which case the continuum limit of 
d^iS) coincides with d{S'^). However, it might also come 
as a string of mainly 4-d cubes at arbitrarily small lattice 
spacing, and then the limit of d^iS) will be significantly 
larger than d{S'^). At the same time, d'"{S) > 1 for the 
connected lattice structures and the physical dimension 
will not be underestimated. 

Since the super-long-distance property involves con- 
nected parts spreading over maximal distances, and since 
it is insensitive to the change of lattice spacing, we can 
conclude that the skeleton is at least 1-dimensional. To 
see whether higher dimensions (less than four) could 
be relevant we consider sets S-^ and their partitions 
= iS-^+ U S-^^ into coherent subsets containing sites 
with positive and negative sign of q^- In Fig. ^ we show 
the /-dependence of the average lattice coherence dimen- 
sion (d^lS-^"^)) for ensembles £i (I) and £2 (11). We see 
that the dimension slightly decreases as we approach the 
continuum limit. The naive linear extrapolation for the 
skeleton value f = 0.18 (pentagon) gives {d-'^{Sf^)) « 1.5 
in the continuum limit and we thus conclude that the 
physical local dimension associated with the skeleton is 
^ ^ d < 1.5. Assuming that the local structure is not 
fractal (an option to be investigated) we are thus led to 
the conclusion that the skeleton is locally 1-dimensional. 

For comparison, we show in Fig.^jalso (d-'^^S-^^)) for 
£2 after random permutation of sites (III), and the re- 
sult for effective density qi^~'^^'^ MoV) same lattice 
spacing (IV). As expected, the curve (III) lies lower and 
tends to zero smoothly at small fractions. On the other 
hand, for the effective density one naively expects 4-d 
coherence on the scale 1/A in the continuum limit, and 
thus the lattice coherence dimension is much higher. 

5. Finally, to start developing a visual intuition for 
the kind of structure we are observing, we now proceed 
to graphically verify the existence of the coherent sheets 
described in part 2. To do that, we have studied the 
graphs of on 2-d sections of the underlying space-time 



torus. The generic behavior of qx is shown in Fig. 5 (con- 
figuration from ensemble £1). One can clearly observe the 
enhancement of the structure with linear "ridges" spread- 
ing over the whole system. It is easy to see that this is 
precisely the behavior expected from the (almost) space- 
filling super-long-distance sheets. Indeed, recall that if 
51 is an arbitrary base manifold and Qi,D,2 C fi, then 
dim(riinf22) = dim(rii)-|-dim(il2) — dim(il) generically. 
If we identify f2 with the 4-d space-time torus, fli with 
the 2-d torus of the section plane, and ^2 with the man- 
ifold of the structure, we get 

dim(r2i n ^2) = dim(r22) - 2 (2) 

Thus if O2 is a 3-d super-long-distance hypersurface, we 
should generically see very extended 1-d regions of sign- 
coherence, which is precisely what is observed. Thus, at 
the lattice cutoff of about 2 GeV the sheets behave as 3-d 
hypersurfaces. The precise local dimension of maximal 
coherent regions in the continuum limit will be studied 
elsewhere. 

In this work we have addressed the following question. 
Can we identify a fundamental (incorporating all scales) 
ordered structure in topological charge density for typical 
configurations contributing to the pure-glue QCD path 
integral? In other words, using nothing more than a local 
operator q^ to calculate space-time distribution of topo- 
logical charge, can one detect well-defined patterns in 
these distributions? While such structure has not been 
observed before, we have argued that the special prop- 
erties of the newly available topological charge densities 
associated with Ginsparg- Wilson fermions give reasons 
to re-examine this question in detail. In particular, the 
suppression of ultraviolet noise by chiral smoothing sug- 
gests that the problem of ultraviolet dominance might 
be solved and the underlying order could be revealed. 
Performing numerical experiments, we have then pro- 
vided evidence for a fundamental structure of this type by 
showing that topological charge exhibits long-range order 
(sign-coherence) on lower-dimensional subsets of the 4-d 
Euclidean space. In this initial study we have concen- 
trated on certain geometrical properties of this structure. 
Indeed, if topological charge fluctuations play a role in 
important aspects of the QCD vacuum (such as SxSB), 
then some geometrical features should reflect that. While 
it is not clear yet which structural aspects are physically 
most relevant, we emphasize two that we find most in- 
triguing: (i) Space-time regions with intense topological 
charge density are organized into a maximally extended 
(super- long-distance) structure ( "skeleton") . This sug- 
gests that there is no specific long-distance scale associ- 
ated with individual geometrical objects. Rather, there 
is a global low-dimensional structure that might have 
to be considered as a whole. (Distance scales associ- 
ated with Aqcd could be reflected by local features of 
such global structure.) We speculate that the super-long- 
distance character of the skeleton might be naturally as- 
sociated with the long-range propagation of Goldstone 
pions. It will be demonstrated elsewhere that the super- 
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long-distance property is inherited and enhanced at low 
energy, (ii) Our data indicates that the skeleton is a 
locally-linear folded structure, much like a protein or a 
neural network in 3-d. One aspect that might be signif- 
icant about this result is that the propagation of mass- 



less quarks on a 1-d manifold naturally leads to a non- 
zero density of Dirac eigenvalues around zero, and hence 
SxSB (see e.g. Ref. 19] for discussion of similar issues). 
The possibility that such a scenario might be relevant for 
QCD is being studied. 
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